Two particle excited states entanglement entropy in a one-dimensional ring 
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The properties of tlie entanglement entropy (EE) of two particle excited states in a one- 
dimensional ring are studied. For a clean system we show analytically that as long as the momenta 
of the two particles are not close, the EE is twice the value of the EE of any single particle state. 
For almost identical momenta the EE is lower than this value. The introduction of disorder is 
numerically shown to lead to a decrease in the median EE of a two particle excited state, while 
interactions (which have no effect for the clean case) mitigate the decrease. For a ring which is of 
the same size as the localization length, interaction increase the EE of a typical two particle excited 
state above the clean system EE value. 

PACS numbers: 73.20.Fz,03.65.Ud,71.10.Pm,73.21.Hb 



INTRODUCTION 

There has been a growing interest in the behavior of 
entanglement entropy (EE) [l| in different physical fields. 
In condensed mater, much of the interest stems from the 
behavior of the ground state EE in the presence of quan- 
tum phase transitions (QPTs) [H-Q. The EE of a finite 
region A of a one-dimensional system grows logarithmi- 
cally as long as the region's size La is smaller than the 
correlation length ^ characterizing the system, while it 
saturates for La > C 043- This behavior may be used 
in order to extract ^, for example the ground state local- 
ization length of the Anderson transition 8]. 

A natural question is what is the behavior of the EE for 
the excited states 0] ? Beyond the growing interest in EE 
coming from the quantum information circles, the ques- 
tion whether EE is a useful concept in studying the be- 
havior of excited states, is relevant to the condensed mat- 
ter community. Low lying excited states in the vicinity of 
a ground state quantum critical point (QCP) should be 
strongly influenced by the critical point [10| , and one ex- 
pects it to show in the behavior of the EE of these states. 
Moreover, the whole concept of the many-body localiza- 
tion transition is centered on the behavior of the 
excited states. The localization-delocalization transition 
occurring at a critical excitation energy should change 
the properties of excitations above it which should be 
manifested in the properties of the excited states. Al- 
though much effort went into trying to understand the 
transition using different properties of the excited states 
(such as level statistics, inverse participation ratio, con- 



ductance, and correlations) (1J-|20|, all these studies were 



performed for rather small systems, and many questions 
remain open. Recently 21| , time evolution of the entan- 
glement of an initial state was studied, and it showed 
signs of many-particle delocalization. Thus, EE seems as 
a useful tool to study the many-body localization transi- 
tion. 

Unlike the ground state EE for which universal results 
exist, the understanding of EE for the excited states is 
still a work in progress 22-23]. Therefore, it would be 



useful to consider a system for which the EE of the ex- 
cited states is simple enough to describe analytically, al- 
though it exhibits interesting behavior such as interaction 
induced delocalization of the excited states. In this paper 
we study the EE for such a system, namely two particles 
on a ring. The study of two interacting particle (TIP) in 
a disordered one dimensional system has a long history 
in the context of many particle delocalization problem. 
All single-electron states for any amount of disorder are 
localized [25j. This continues to be true for two-electron 
states, however, the localization length becomes longer as 
the repulsive interaction becomes stronger [2^ 27 1 . This 



interaction induced delocalization was confirmed numer- 



ically |28l - l3l( . It is important to emphasize that there is 
no enhancement of the localization length in the ground 
state. The delocalization becomes significant only for 
higher excitations. 



CLEAN RING 

For a clean ring composed of iV sites, the tight-binding 
Hamiltonian is given by: 



H 



N 



N 



t^(e^"a]a,+i-|-/i.c.), (1) 
i=i 



where for the clean case Cj = 0, i = 1 is the hopping 
matrix element between neighboring sites, and Oj is the 
creation operator of a spinless electron at site j on the 
ring. In order to break symmetry (we shall see why this 
is important further on) a magnetic flux threading the 
ring is introduced, where a — 27r0/ ((/)qN), and 4>o = hc/e 
is the quantum flux unit. The single-electron eigenvalues 
are e{k) = —2tcos{p — a) where p = 2Trk/N, and k — 
0,±1,±2,...± N/2. The eigenvectors are given by 



N 



|fc) = (l/%/7V)^exp(zpj>J|0), 



(2) 



where 10) is the vacuum state. 
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The two particle eigenvalues are e{kl,k2) = 
—2t (cos(pi — a) + cos(p2 — a)). The eigenvector 

N 

|fei,fe) = (l/iV) A{k,,n,k2,j2)a]^a]^\%), (3) 

ii>j2=i 

where 

A{ki,ji,k2,32) = 

{exp(i{piji +P2j2)) - exp{i{pij2+P2ji)))- (4) 

Once the eigenvectors of the system are available one 
can (in principal) calculate the EE. The entanglement 
between a region A (of length Na) in the system and the 
rest of the system (denoted by B) for a given eigenstate 
1 5') is measured by the EE Sa/b- This EE is related to 
the region's reduced density matrix pa/B: defined in the 
following way: 

PA/B = TrB/A\'f){-9l (5) 

where the trace is over region's B or A degrees of freedom. 
The EE is related to the eigenvalues of the reduced 
density matrix: 

Sa/b = -S^Ai In(Ai). (6) 

One important result of this definition is the symmetry 
between the EE of the two regions Sa ~ Sb- 

Following Cheong and Henley [s^l, one can write the 
pure state |^) — ^ \'iA)\(l>B,i), where is a complete 

i 

orthonormal many-body basis of region A, while \4>B,i) 
is the state in region B associated with Please note 
that \(f>B,i) is not normalized. Using that notation, the 
reduced density matrix, 

N 

PA^ \^A){jA\, (7) 

or in matrix form: 

PA{i,j) ^ {(I^bMbj)- (8) 

Utilizing the occupation basis in the A region, i.e., \iA) = 
\n\,n2,nl, . . . , n]^^) (where = 0, 1), one can define an 

Na 

operator, Ki = + (1 ^ "i^D^sO-l], resulting in: 

PAit,j)^{m}Kj\^)- (9) 

It is important to note that pA{i,j) ^ only for states 
which have the same number of particles n\ in region 

N 

A, where n\ — n\. Thus, in this basis, the reduced 

s=l 

density matrix is composed of blocks which increases in 



size with n\. Thus for n\ = 0, the block size is one, for 
= 1 it is Na, for = 2 it is Na{Na - l)/2, etc. 

Thus, the task of calculating the EE of a region A at a 
given excitation 1 5') is equivalent to calculating the eigen- 
values of the matrix pA(i,j)- Since the blocks are uncou- 
pled, it is possible to diagonalize each block with a given 
number of particles, p^''(i, j), independently (where £ de- 
notes the number of particles in region A). For a state 
1 5*), which is the ground state of a half- filled ring at 
N oo, it is possible (in several different ways) to show 
that Sa = -(l/3)ln(a;) + Const Q, where x = Na/N. 
For the excited states the task becomes more difficult, 
and no simple and general result for Sa exists Q. Here 
we will calculate the Sa for two-particle excitations. 

For the sake of completeness lets first consider single 
particle state entanglement. In this simple case pA{i, j) is 
composed of two blocks: p^^(l,l) and p^^{i,j) (where 
i,j = 1...Na)- Direct evaluation of Eq. ^ for any 
single-particle state = \k) results in p^''(l, \) — 1 — x, 
while p^^(i,j) = (1/7V) exp(— ip(i — j)). The latter is 
a Toeplitz matrix, with one eigenvalue equal to x and 
Na — 1 zero eigenvalues. As expected, the EE for any 
single particle eigenstate |fc) is equal to 

S'a = -a;ln(a;) - (1 -a;)ln(l -x) (10) 

and does not depend on |fc). 

A note of caution is in place. Since the eigenvalues 
for k and —k are degenerate, any linear combination of 
|fc) and I — k) are an excited state of Eq. ([!]). These 
linear combinations have different values of the EE, and 
therefore, strictly speaking, the EE for degenerate ex- 
cited states is ill defined. We circumvent this problem by 
introducing a degeneracy breaking magnetic flux (j) into 
Eq. ([1]). As long as the degeneracy is broken the EE 
of any excited state ^2) is well defined and docs not 
depend on (j). 

For two-particle states |fci,fc2), the reduced density 
matrix is composed of three blocks: /9^'(1,1), PA\hj) 
(of size Na) and /9^(i,j) (size Na{Na - l)/2). For the 
zero particle block: 

— Y \A{ji,k,,j2,k2)\'^{l-xr-y', (11) 
ji>n>NA 

where y = (sin(7r(fc2 — ki)x)/ {TT{k2 — ki))). Thus the 
eigenvalue of this block is (1 — x)^ — y^. Using symme- 
try, one can immediately deduce the eigenvalues of the 
two-particle reduced density matrix, without actually di- 

agonalizing the 7Vyi(iV^ — 1)/2 matrix. Since Sa = Sb the 

(2) 

contribution to the EE from (z,j), must be equal to 
the contribution from p^-*(l, 1). This infers that p^2\hj) 
has only one non-zero eigenvalue. Seeing that region B's 



length is A'' — Na, according to Eq. (|lip . the non-zero 
eigenvahie of p^'^{i,j) is equal to — y^. 

The one-particle block density matrix is given by: 



1 ^ 

(*'-^) = lv2 X! ^*(il''=l'*'^2)^(jl,fcl,j, ^2) 



iV2 
-jpi(i-i) 

1 



lN{p2 - Pi 



]i>Na 



. |^g-*(P2-Pl)ij-g«(P2-Pl)JVA _ 
,»(P2-pi)j(-g-i(p2-pi)AfA — 1)]^. (12) 



This cumbersome form is substantially simplified when 
k2 — ki is large. In that case the second term in Eq. 
(fT2)l may be neglected and the density matrix block has 
a Toeplitz form 



(13) 



with Na — 2 zero eigenvalues, and two degenerate eigen- 
values equal to x{l — x). The second term is negligible 
also when x — Na/N ~ 1/2, and ^2 — fci is even, result- 
ing in the same eigenvalues. We do not have the gen- 
eral solution, nevertheless, it can be shown numerically 
that p^'(j, j) has no more than two non-zero eigenvalues, 
which depend only on the difference k2 — ki. Moreover, 
since the sum of all eigenvectors of the density matrix 
should be one, the sum of those two eigenvalues should 
be 2a;(l — x) -I- 2y^. For the ground state (and excita- 
tions for which ^2 — fci = 1) the two eigenvalues are well 
described by 2a;(l — a;) -I- (2 — l/n)y'^ and j/^/tt. 

Thus, the EE of two-particle states composed from two 
single-particle states of significantly different wave num- 
bers which are the majority of the two-particle states. 



SA{k2 - fci > 1) = -2[(1 -x) ln(l - a;) -xln(a;)] .(14) 



This is twice the EE of a single particle state (Eq. (ITOl) . 
Thus, as long as the two occupied states |fci) and 1^2) 
are far enough from each other, the two-particle EE is 
just the sum of the EE of each occupied state. In the 
opposite limit 

SA{k2 - fci = 1) ^ -((1 - xf - 2/2) ln((l - x)2 - y2) 
- (x2 - ln(.x2 _ y2) _ [^2/^] i^[y2/^] 

- [2x(l - x) + (2 - l/7r)y2] ln[2x(l - x) + (2 - l/n)y% 

(15) 

The EE curves for other values of ^2 — fci can be calcu- 
lated numerically by diagonalizing the Na x Na matrix 
representing p^'' (the two other eigenvalues for and 
Pj^ are given be Eq. ([TT]) '). The results are depicted in 




FIG. 1: (Color online) The EE of a two particle state |fci, ^2) 
of a clean system as function of region A's size x = Na/N for 
a system of length A'' = 1000. The Heavy lines correspond 
to the analytic prediction, (dotted for k2 — ki = 1, Eq. (|15|) . 
dashed for k2 — k\ 3> 1, Eq. HU). The thin lines pertain to 
the numerically calculated EE for all values of ^2 — k\ between 
1 . . . 30, where odd values are depicted by red curves, and even 
one by red curves. It is clear that for k'z — fei larger than 5 the 
numerical curves fit Eq. (|14|) quite well. It is also clear that 
all even k2 — k\ reach the same EE value Sa{x = 1/2) = ln(4) 
once Na = N/2. 



Fig. [TJ All the values of the EE for any value of k2 — ki 
fall in between those two limits, but as can be seen in Fig. 
[T]they quite quickly fall on the A:2 ~ ^1 ^ 1 curve. Since 
there is a large phase space for fc2 — fci ^ 1, a typical 
two particle excitation corresponds to the EE described 
in Eq. [H 

Another interesting behavior that can be gleaned from 
Fig. [l]is that all two-particle states of even k2 — ki reach 
the same EE value at Na — N/2. This stems from the 
structure of p^^(z, j) (Eq. (fT^ l. where the two last terms 
are multiplied by e'^^p^^p^')'^^ - 1 = e-"('=2-'=i) -1 = 
(since P2—P1 = 2n{k2 — ki)/N) returning to the density 
matrix block with a Toeplitz form depicted in Eq. (|13l) . 
and the corresponding two degenerate eigenvalues equal 
to 1/4. Since at Na = N/2, x = 1/2, and y = 0, the two 
other block eigenvalues are also 1/4, resulting in SA{k2 — 
ki = odd,x = 1/2) = ln(4). Thus the largest EE in a 
two-particle clean ring system is equal to ln(4). 



INTERACTING CLEAN RING 

Incorporating nearest neighbor electron-electron inter- 
actions into the system, results in adding an interaction 
term given by 



N 



(16) 



to the Hamiltonian, iJ, depicted in Eq. ([T]). In a clean 
system it is well known that far from half-filling the sys - 
tem behaves as a Luttinger liquid for any value of U [3J] . 
For the ground state EE of a clean system at half-filling 
(and U <2, i.e., a Luttinger liquid) the EE changes only 
by an overall constant [l[ Q, while retaining the same 
logarithmic dependence. Thus, we expect that the EE 
of the two-particle states in a clean system will not be 
essentially affected by the presence or absence of electron- 
electron interactions. Unfortunately, it is not possible to 
calculate analytically the two-particle states of the in- 
teracting system. Thus, we must rely on a numerical 
solution for the problem. 

Exact diagonalization is used to calculate all the 
eigevectors of H + Hint, represented by a N{N — l)/2 
matrix. We have chosen a 100 site system, resulting in 
a matrix of size 4950. A reduced density matrix pA, of 
size 1 -I- Na + Na{Na + l)/2 is then constructed and di- 
agonalized for each eigenstate, and the EE is calculated 
using its eigenvectors according to Eq. The results 
are shown in Fig. ^ where the EE of 31 states around 
the ground state (i.e. the ground state and 1st - 30th ex- 
citation), and at quarter of the two-particle band (1222th 
- 1252st excitation) are shown. In both cases the EE for 
the non-interacting {U = 0) as well as for the interacting 
{U = 1) cases are almost equal (the interacting case is 
larger by a minute constant (of order lO^'' which can not 
be resolved at the resolution of the figure). As expected 
around the ground-state the excitations belong to the low 
^2 — fci sector while for the higher excitations most states 
corresponds to large values of k2 — ki, i.e., well described 
byEq. dm). 



INTERACTING DISORDERED RING 

When disorder is added to a non-interacting system, 
all single particle states become localized. For the many- 
particle states, the behavior is more involved. As long 
as no interaction is present, the many-particle states re- 
main localized, both for the ground state [s^ as well as 
for all the excited states. Once interaction is introduced, 
the ground-state as well as low lying excitations remain 
localized, while above a critical energy the many-particle 
excitations are predicted to delocalize 11- 1!]. This tran- 



sition, termed the many-body or Fock space localization 
transition, stems from the interactions coupling excita- 
tions with a different number of electron-hole couples. 
This type of transition is irrelevant for two particle sys- 
tems. Nevertheless, as argued by Shepelynsky and Imry 
2^ 27 1, interaction between the pair of particles should 
enhance the two particle localization length, compared 
to the single electron localization length, as long as the 
two-particle level spacing is significantly smaller than the 
single electron level spacing, i.e., for higher excitations. 
Can we see any signature of the enhanced two-particle 
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FIG. 2: (Color online) The EE of 31 states of a clean system as 
function of region A's size Na for a system of length — ICQ, 
in the absence {U — 0, black line), or presence {U — 1, dotted 
red line) of electron-electron interactions. Panel (a) depicts 
states in the vicinity of the ground-state (the ground state and 
1st - 30th excitation), panel (b) shows the excitations around 
quarter of the two-particle band (1222th - 1252st excitation). 
The Heavy dashed line correspond to the analytic prediction 
for ^2 — fci ^ 1, Eq. HU). For a clean system the interaction 
has no influence on the EE. 



localization length in the EE behavior of the excited 
states? First we have to understand the influence of dis- 
order on the EE. In Ref. [^it has been shown that for the 
ground-state the EE saturates on the length scale of ^, 
and does not continue to grow logarithmic as in a clean 
system. Thus, the EE of a disordered system is always 
lower than the EE of a clean system. One would expect 
this feature to hold also for excited states. We check this 
assumption by calculating the EE using the excitations 
of the Hamiltonian given in Eq. ([1]), where the disorder 
is represented by a random on-site energy, ej taken from 
a uniform distribution in the range [—W/2,W/2]. For 
W = 3, single electron states at the middle of the band 
are expected to have a localization length ^ ~ 10 [sHi, 
while close to the band edge the single electron states 
are supposed to be much more strongly localized (Lif- 
shitz tails) (s^j- 

The EE is calculated by exact diagonalization for sys- 
tems of size N = 100 as described in the previous section. 
The results are presented in Fig. [3J where the median EE 
in the vicinity of the ground-state (1st - 30th excitation), 
at 1/16 of the band (294th - 324th excitation), at 1/8 of 
the band (603th - 633th excitation), and at 1/4 of the 
band (1222th - 1252th excitation). The median EE is 
taken across the 31 excitations in each segment and 50 
different realizations of disorder. We calculate the EE in 
the absence [U = 0) and presence {U = 1, U = 3) of 
electron-electron interactions. Around the ground state, 
the interactions do not play a significant role, and for all 
cases the EE is strongly suppressed compared to typical 



1.4 



1.2 



-<0.8 



CO 



' I ' 

G — o arouiig GS U=0 
G-^ around GS U=l 
around GS U=3 

□ ■ □ 1/16 band U=0 

□ ■ □ 1/16 band U=l 
D □ 1/16 band U=3 
«-« 1/8 band U=0 
0--« 1/8 band U=l 

- 1/8 band U=3 
i -A 1/4 band U=0 
i -A 1/4 band U=l 

- 1/4 band U=3 
-- Clean 




FIG. 3: (Color online) The median EE of an ensemble of 
states for different regions of the band collected from 20 dif- 
ferent realizations of strong disordered {W = 3), as function 
of region A's size Na for a system of length A'^ = 100. Contin- 
uous lines depict the median EE in the vicinity of the ground- 
state (1st - 30th excitation), dotted lines at al/16 of the band 
(294th - 324th excitation), dashed lines at al/8 of the band 
(603th - 633th excitation), and dot-dashed lines at al/4 of the 
band (1222th - 1252th excitation). Black lines correspond to 
no electron-electron interactions, while red lines indicate the 
presence of moderate electron-electron interactions {U = 1), 
and green lines correspond to stronger interactions {U = 3). 
The heavy dashed line correspond to the maximum EE for a 
clean system (Eq. (|14|l 'l. Error bar represent a range between 
the 40th and 60th percentiles. 



values for a clean system. This is expected, since as men- 
tioned, at the bottom of the band the states are strongly 
localized, and therefore the EE is low. 

For the non-interacting case, the EE higher in the band 
are less suppressed as the localization length grows. The 
EE for high excitations in the interacting case is always 
larger than for the corresponding non-interacting states. 
This is a clear signature for the effect of interactions on lo- 
calization of the two-particle states, which become more 
entangled as interaction is present, although there is no 
significant difference between U = 1 and U = 3. It is 
also clear that for higher excitations (larger localization 
length) the enhancement of the EE becomes stronger. 

This enhancement could b e expe cted on physical 
grounds. 
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FIG. 4: (Color online) As in Fig. Ofor weak disordered {W = 
1). For high enough excitation energy, the median EE in the 
presence of disorder and interactions exceeds the maximum 
clean value indicated by the heavy dashed line (Eq. p4p '). 



der. Surprisingly, above 1/8 of the band (corresponding 
to an excitation energy of t), the EE of the disordered 
interacting system is significantly larger than the limit 
for a clean system (ln(4)). Although, extrapolating from 
the results presented in Fig. [31 increasing the system size 
while keeping the disorder fixed will result in a decrease 
of the EE below the clean system values once L ^ ^. The 
increase above the clean system excitation EE may stem 
from the fact that as long as the two particles are con- 
fined within a single localization length the particles can 
not avoid each other and spend much time close to each 
other, leading to an enhancement of the EE. When the 
system size is much larger than the localization length, 
the two particles can reside in different regions of the 
sample, and interactions will not play an important role. 
However, this hand waving picture requires further study. 

At first glance, these results seem to indicate that 
although interactions may enhance the EE as long as 
L < ^, they become irrelevant for L ^ ^, showin g no 



261- 31|, the localization 



As has been shown 
length associated with an interacting two-particle state 
is larger than for a non-interacting state with the same 
disorder. Thus, one expects that the EE will also be 
larger and closer to its clean system value. 

We therefore also investigate the case of weaker disor- 
der, for which the localization length is of order of the sys- 
tem size {W = 1, ^ ^ 100). As can be seen in Fig. HI for 
the non-interacting case a similar pattern to the one ob- 
served Fig. [3] remains, although the EE is less suppressed 
by the weaker disorder. As expected, the enhancement 
of the EE by interactions is stronger for the weaker disor- theless, the fact that the two-particle behavior fits nicely 



sup port for the many-particle delocalization scenario [11 
|l3| which should occur for L ^ ^. This interpretation 
is wrong, since the many-particle delocalization scenario 
deals with a constant density of particles, and delocal- 
ization is predicted only when there are at least a couple 
of particles in the range of a single particle localization 
length. Thus the observed two-particle EE enhancement 
when the two particles are within a distance of ^, as well 
as the fact that the enhancement increases significantly 
when the excitation energy, fits nicely with the scenario 
promoted in Ref. |l3|. Of course, coupling between states 
with a different number of electron- hole generation is cru- 
cial for the delocalization scenario, and therefore a full 
demonstration of the delocalization transition has to be 
performed for a finite electron density system. Never- 



6 



with the delocahzation scenario is encouraging. 



CONCLUSIONS 

The properties of the EE of two particle excited states 

in a one-dimensional ring were studied. For a clean sys- 
tem, the EE depends only on the difference in momentum 
between the two particles. If the difference is large the 
EE corresponds to the EE of two independent single par- 
ticle states, i.e., Sa = — 2[a;ln(a;) + {1 — x) ln(l — x)]. On 
the other hand, if the momenta are close, the EE of the 
two particle state is reduced compared to this value. 

One may extrapolate that for m particles on a iV site 
ring, as long as the density is low {m/N ^1), the upper 
limit of the EE is Sa = — m[a;ln(a;) + {1 — x) ln(l — x)], 
which is also the typical value. This will be valid if the 
difference between the momenta of all particles taking 
part in a particular many-particle excited state is large. 
If this is not the case we expect the EE of the excited 
state to be lower. Further investigation of these cases is 
underway. 

We have verified numerically that disorder reduces the 
EE. Short range particle interaction leads to an enhance- 
ment of the excited state EE, which become very signifi- 
cant once the localization length is of order of the system 
size. For high excitations the median EE of a many parti- 
cle interacting excitation is not only above the disordered 
case, but exceeds the clean system limit. This may be 
related to the fact that localization forces the two par- 
ticles to dynamically spend more time in the vicinity of 
each other, although this argument merits further study. 
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(Grant 686/10) is gratefully acknowledged. 



[1] For recent reviews see: L. Amico, R. Fazio, A. Osterloh, 
and V. Vedral, Rev. Mod. Phys. 80, 517 (2008); J. Eisert, 
M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 277 
(2010); and references therein. 

[2] M. Vojta, Pliil. Mag. 86, 1807 (2006). 

[3] K. Le Hur, Ann. Phys. 323, 2208 (2008). 

[4] M. Goldstein, Y. Gefen and R. Berkovits, Phys. Rev. B 
83, 245112 (2011). 

[5] C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. B 
424, 443 (1994). 

[6] G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. 
Rev. Lett. 90, 227902 (2003); J. I. Latorre, E. Rico, and 
G. Vidal, Quant. Inf. Comp. 4, 048 (2004). 

[7] P. Calabrese and J. Cardy, J. Stat. Mech. P06002 (2004). 



[9 

[lo; 
[11 

[12 
[13 
[14 

[15 

[16 

[17 

[18 
[19 

[20 

[21 

[22; 

[23 
[24; 

[25; 

[26' 
[27 
[28; 

[29; 

[3o; 

[31 

[32 

[33 
[34 

[35; 

[36 
[37 



R. Berkovits, Phys. Rev. Lett. 108, 176803 (2012). 

V. Alba, M. Fagotti and P. Calabrese, J. Stat. Mech. 

P10020 (2009). 

S. Sachdev Quantum Phase Transitions Cambridge Uni- 
versity Press (1999). 

B. L. Altshuler, Y. Gefen, A. Kamcnev, and L. S. Levitov, 
Phys. Rev. Lett. 78, 2803 (1997). 

I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov Phys. 
Rev. Lett. 95, 206603 (2005). 

D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Ann. 
Phys. (N.Y.) 321, 1126 (2006). 

R. Berkovits and Y. Avishai, Phys. Rev. Lett. 80, 568 
(1998); R. Berkovits and B. I. Shklovskii, J. Phys. Con- 
dons. Matter 11, 779 (1999). 

V. Oganesyan and D. A. Huse, Phys. Rev. B 75, 155111 

(2007) ; V. Oganesyan, A. Pal, and D. A. Huse, ibid. 80, 
115104 (2009). 

C. Monthus and T. Garel, Phys. Rev. B 81, 134202 

(2010). 

T. C. Bcrkclbach and D. R. Reichman, Phys. Rev. B 81, 

224429 (2010). 

A. Pal and D. A. Huse, Phys. Rev. B 82, 174411 (2010). 

E. Canovi, D. Rossini, R. Fazio, G. E. Santoro, and A. 
Silva, Phys. Rev. B 83, 094431 (2011). 

E. Cuevas, M. Fcigel'man, L. loffe and M. Mezard. Na^ 
ture Communications 3, 1128 (2012). 

J. H. Bardarson, F. PoUmann, and J. E. Moore, Phys. 
Rev. Lett. 109, 017202 (2012). 

F. C. Alcaraz and M. S. Sarandy, Phys. Rev. A 78 032319 

(2008) . 

L. Masanos, Phys. Rev. A 80 052104 (2009). 

M. I. Bcrganza, F. C. Alcaraz, and G. Sierra, J. Stat. 

Mech. P01016 (2012). 

For a review sec: P. A. Lee and T. V. Ramakrishnan, 
Rev. Mod. Phys. 57, 287 (1985). 

D. L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994). 
Y. Imry, Europhys. Lett. 30, 405 (1995). 

K. Frahm, A. Miillcr-Grocling, J.- L. Pichard, and D. 

Weinmann, Europhys. Lett. 31, 169 (1995). 

D. Weinmann, A. Miiller-Groeling, J.- L. Pichard, and 

K. Frahm, Phys. Rev. Lett. 75, 1598 (1995). 

F. von Oppen, T. Wettig, and J. Miiller, Phys. Rev. Lett. 

76, 491 (1996). 

Ph. Jacquod, D. L. Shepelyansky, and O. P. Sushkov, 
Phys. Rev. Lett. 78, 923 (1997). 

S.-A. Chcong and C. L. Henley, Physical Review B, 69, 
075111 (2004). 

L. Gong and P. Tong, Phys. Rev. B 76, 085121 (2007). 

T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, New York, 2003). 
W. Apcl, J. Phys. C 15, 1973 (1982); W. Apcl and T. M. 
Rice, Phys. Rev. B 26, 7063 (1982); T. Giamarchi and H. 
J. Schulz, Phys. Rev. B 37, 325 (1988). 
R. A. Romer and M. Schreiber, Phys. Rev. Lett. 78, 515 
(1997). 

T. M. Lifshitz, Adv. Phys. 13, 483 (1964). 



